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We present a simple extension of the standard model that gives rise to baryogenesis a has a dark 
matter candidate of ©(GeV) mass. A minimal set of new fields required for baryogenesis includes 
two ©(TeV) colored scalars and a singlet fermion. The fermion also becomes a viable dark matter 
candidate when it is nearly degenerate in mass with the proton. Dark matter and baryon asymmetry 
are produced form the decay of heavy scalars, which can lead to a natural explanation of the baryon- 
dark matter coincidence problem. The dark matter candidate escapes direct and indirect detection, 
but can be probed at the LHC. The supersymmetric extension of this model is straightforward and 
leads to a multi-component dark matter scenario, which improves the direct and indirect detection 
prospects. 



I. INTRODUCTION 

Weakly interacting massive particles (WIMPs) are 
promising dark matter (DM) candidates [l|. WIMPs 
typically arise in models of particle physics beyond the 
standard model (SM). In supersymmetric (SUSY) mod- 
els with conserved -R-parity, the lightest supersymmetric 
particle (LSP) is the DM candidate. In the standard 
scenario, the DM relic abundance, precisely measured 
by cosmic microwave background experiments |^, is ex- 
plained via thermal freeze-out of LSP annihilation in the 
early universe. 

However, the current LHC bounds (3| have put this 
scenario under increasing pressure. The fact that no 
SUSY particles have been found so far, keeps pushing 
up the limits on the mass of gluinos and squarks of the 
first two generations, and to a lesser extent the stop and 
sbottom. This has motivated new scenarios, like natural 
SUSY [4l|, that can accommodate heavy SUSY particles 
in accordance with the Higgs mass. Although the masses 
of the non-colored particles do not have much of a con- 
straint, but based on the allowed parameter space of the 
new scenarios, thermal DM is not favored unless the DM 
particle is very heavy too. For example, within natural 
SUSY scenarios the Higgsino typically arises as the DM 
candidate. The annihilation rate in this case is larger 
than the nominal value for thermal scenario 3 x 10 ~^^ 
cvD? /sec for sub-TeV DM mass. On the other hand, such 
large annihilation rates are becoming more and more con- 
strained by the Fermi-LAT data from DM annihilation in 
the dwarf galaxies and the galactic center 'S]. As a result, 
one requires non-thermal mechanisms and/or other DM 
candidates in order to obtain the correct DM abundance. 

Another challenge for physics beyond the SM is the 
explanation of the baryon asymmetry of the universe [6| . 
Constraints on the stop mass have also put the elec- 
troweak baryogenesis [7| in the minimal supersymmetric 
standard model (MSSM) in a tight corner. There are 
many alternative scenarios that can explain the matter- 
antimatter asymmetry like leptogenesis Q, Affleck-Dine 



baryogenesis [0], hidden sector baryogenesis [13], etc. 
One curious observation is that the energy densities in 
baryons and DM are of the same order of magnitude, the 
so-called "baryon-DM coincidence puzzle" . Then a ques- 
tion arises as whether this apparent coincidence may be 
addressed by an underlying connection between the DM 
production and baryogenesis scenarios |ll| . 

In this work, we present a minimal extension of the SM 
to address these questions. We introduce renormalizable 
baryon number violating interactions in the Lagrangian 
that can lead to a successful baryogenesis. The mini- 
mal field content that is required to achieve this includes 
iso-singlet color-triplet scalars and one singlet fermion. 
We show that the fermion becomes stable, hence a DM 
candidate, when its mass is around O(GeV). The DM 
relic density and the baryon asymmetry are produced 
non-thermally from the decay of some heavy particle(s). 
Non-thermal baryogenesis has the virtue that couplings 
associated with the new fields do not need to be artifi- 
cially small. Moreover, the non-thermal mechanism can 
correlate the DM relic abundance and baryon asymme- 
try. Since DM mass is C'(GeV), a correlation between 
the number densities will automatically translates into 
a similar relation between the DM and baryon energy 
densities. This can provide a natural explanation of the 
baryon-DM coincidence puzzle. 

In this model, the DM candidate interacts with up- 
type quarks via the exchange of colored scalar fields. 
We see that the resulting spin-independent and spin- 
dependent DM-nucleon scattering cross sections are well 
below the bounds from current and upcoming experi- 
ments, which makes the prospects for direct detection 
weak. In addition, due to its low mass of ©(GeV), in- 
direct signals from DM annihilation will be negligible. 
However, the model may be probed at the LHC via the 
colored scalars if they have ©(TeV) masses. The de- 
cay channel including the DM candidate will give rise 
to a missing energy signal. The SUSY extension of this 
model is straightforward. In this case the scalar partner 
of the DM can also become stable, if i?-parity is con- 



served, which allows a scenario with multi-component 
DM. 

This paper is organized as follows. In Section II, we 
present the model and see how it can lead to a DM candi- 
date that is nearly degenerate with the proton. In Section 

III, we discuss non-thermal production of DM and baryo- 
genesis in the model, and show how this can naturally 
address the baryon-DM coincidence problem. In Section 

IV, we discuss prospects for probing the model via direct 
detection experiments and at the LHC, and comment on 
the SUSY extension of the model and the possibility of 
having a multi-component DM scenario. We conclude 
the paper in Section V. 



II. THE MODEL 

We start with the SM Lagrangian and add renormal- 
izable terms to it that violate baryon number in order to 
find a successful baryogenesis scenario. Then gauge in- 
variance requires introducing new colored fields. A min- 
imal set up includes two iso-singlet color-triplet scalars 
Xa (a = 1, 2) with hypercharge -1-4/3. This allows us 
to have the baryon number violating interaction terms 
XacFdP in the Lagrangian (we have used two-component 
Weyl fermions). ^ 

We note that at least two X are needed to produce 
a baryon asymmetry from the interference of tree-level 
and one-loop diagrams in a decay process governed by 
the Xad'^d'^ interactions. However, although necessary, 
this is not sufficient. The reason being that the total 
asymmetry vanishes after summing over all flavors of df^ 
in the final and intermediate states [l^l- One therefore 
requires additional baryon number violating interactions, 
and the simplest renormalizable term as such is X*Nu'^ 
where iV is a SM singlet. This leads to the following 
Lagrangian 
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c^ 
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{Xo^.XlNu'i -f A^,X,d^d5 + -itinNN + h.c.) 



m^jXap + (kinetic terms) 



(1) 



Here i, j denote flavor indices (color indices are omitted 
for simplicity). We note that A^- is antisymmetric under 
i O J. We will discuss the generation of baryon asym- 
metry form this lagrangian in detail in the next section. 
The singlet N , which plays an important role in baryo- 
gensis, can have a gauge charge under a higher ranked 
symmetry group that includes the SM. 

Assuming that tun ^ rna, one finds an effective four- 
fermion interaction Nu'id'^d'f, after integrating out Xa. 



^ One can also add a term Xu'^cF to the Lagrangian if X has the 
same gauge charges as dp (i.e., iso-singlet color triplet with hy- 
percharge —2/3). This leads to similar consequences, and hence 
we do not consider it as a separate case here. 



The existence of this term implies that N decays to three 
quarks if niM ^ ©(GeV). Also, the decay modes N — > 
p-\- e~ +Ve, N ^ p + e^ + Ve are open as long as wln > 
nip + TTie, where rup and me are the proton mass and the 
electron mass respectively. 

However, N becomes absolutely stable if ttin < nip — 
rrie- The important point to note is that stability of N 
is not related to any new symmetry. All required is that 
the proton be stable. Then the same symmetry that 
ensures stability of the proton, combined with the kine- 
matic condition ijin < rup — me, will lead to N being 
a stable particle. This is the first remarkable property 
of the DM in this model: the field N that is required 
for baryogenesis also becomes a DM candidate only if its 
mass is ^ ©(GeV). 

In addition, we note that the decays p ^ N + e^ + 
Ve (pe) are kinematically allowed if mp > rriN + nie. This 
is unacceptable as it will result in catastrophic proton 
decay. Therefore a viable scenario arises provided that 



< mpf < nip + nig 



(2) 



This is the second remarkable property of the model: the 
Lagrangian in Eq. ([1]) gives rise to a viable DM candidate 
if and only if rrijv ~ rUp. 

We therefore see that the model in Eq. ([T]) not only 
gives rise to baryogenesis, but can also yield a light DM 
candidate of ©(GeV) mass. As we will see in the next sec- 
tion, this can lead to a natural explanation of the baryon- 
DM coincidence problem. 

Some comments are in order. First, the smallness of 
toat is protected against quantum corrections since N 
is a fermion. The one-loop corrections arising from the 
XX*Nu'' coupling result in Stun ~ (A/47r)^ln(A/mjf), 
where A is a cut off. This implies that even for A ^ 0{l) 
we have Stun <C ttin. Therefore the relation ttin ~ nip, 
once satisfied at the tree-level, will not be destabilized 
by radiative corrections. In fact, for A ^ 0(0.1) we have 
SniN <, me- 

Second, being a SM singlet Weyl fermion, one may be 
tempted to identify TV with the right-handed (RH) neu- 
trino. This, however, will allow Dirac Yukawa couplings 
of the neutrino HNL in the Lagrangian. Together with 
the terms X*Nu'^ and Xd'^d'^, this induces the dimension- 
7 operator HLu'^d'^d'^, which will lead to a rapid proton 
decay. In order to ensure stability of the proton, one 
therefore needs to forbid the HNL term. This may be 
achieved, for example, by introducing a gauged U{1)l 
symmetry. Then, similar to the case with a gauged 
U{1)b-l, anomaly cancellation will require the existence 
of three RH neutrinos. The U{1)l symmetry does not 
allow a coupling between the RH neutrinos, which carry 
lepton number, and the X*m° combination. Therefore, 
the singlet N that participates in the X*Nu'^ interaction 
term will be decoupled from the lepton sector. This guar- 
antees the stability of the proton, which also implies the 
stability of DM in our model. 

Finally, we note that it is possible to make the DM 
mass much lower than the proton mass by forbidding the 



Xcf^df^ term in the Lagrangian ([T|) by invoking some dis- 
crete symmetry. In this case the proton decay constraint 
does not apply. However, the absence of the Xcfcf term 
also implies that there will be no baryogenesis. 



III. BARYOGENESIS AND DARK MATTER 
PRODUCTION 



In this section, we discuss baryogenesis and production 
of DM in the model. We see that a single decay process is 
the non-thermal origin of both and results in comparable 
abundances of the baryon asymmetry and DM. Consid- 



ering the DM mass m^v 



this can naturally address 



the baryon-DM coincidence puzzle. 



A. Generation of the baryon asymmetry 



their annihilation into gluons as the temperature drops 
below mi^2- This implies that nxn nx2 are too sup- 
pressed at the time of decay to yield the desired baryon 
asymmetry -qs ~ 0(10"^°). 

Therefore, a non-thermal scenario is needed for suc- 
cessful realization of baryogenesis in this model. In a 
possible scenario Xi and X2 are produced from the late 
decay of a scalar field S with mass mg that reheats the 
universe to a low temperature Tr [13l |. Then the overall 
baryon asymmetry will be 



Vb 



Ys ■ — 



1 

8^ 



1 






ml Bri 



-f 



m| Br2 



E.,, IA'h.P + Efe lAifeP E.,, |A'2,,P + Efc IA2.P 



(6) 



The model in Eq. ([T]) can give rise to generation of 
baryon asymmetry through the decay of colored scalars 
Xa- The interference between tree- level and one- loop 
self-energy diagrams will result in the following baryon 
asymmetry per decay of Xi and X2 respectively [13 
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where 



^s{x) 



We note some differences between the asymmetry param- 
eter in this case and that in the leptogenesis scenario [ij] . 
First, there are no one-loop vertex diagrams here. Sec- 
ond, the numerator of the expression for J^{x) ^ con- 
tains X instead of y^. The reason being that the baryon 
number violation in the self-energy diagram arises from 
the couplings X*Nu'^ and Xd'^d'^ in this case, while a 
Majorana mass term is responsible for the asymmetry in 
leptogenesis. 

The generated baryon asymmetry normalized by the 
entropy density s is given by 



riB 



riB 



ei he2 • (0) 

s s 



A natural choice of parameters is |A|, |A'| ^ 0{1) and 
CP violating phases of 0(1), which for mi ^ m,2 results 
in £1,2 ~ 0(0. 1). However, for mi,2 ~ O(TeV), the third 
Sakharov condition [l^ for generating a baryon asym- 
metry (i.e., out-of-equilibrium decay of Xi, X2) requires 
that |A|, |A'| ^ 1 if thermal initial condition for Xi, X2 is 
assumed. Moreover, in this picture, the comoving num- 
ber density of Xi, X2 exponentially decreases due to 



Here Ys = 2>T^/Am.s is the dilution factor due to entropy 
release by the late decaying scalar field S, and Br 1,2 de- 
note the branching ratios for producing Xi and X2 from 
S decay respectively. 



B. Non-thermal production of dark matter 

The DM candidate N reaches equilibrium with the pri- 
mordial plasma at sufficiently high temperatures through 
its interactions with u'^ given in Eq. ([T]). At tempera- 
tures T <C ?Tii,2, the cross section for N scattering off 
u"^, or its pair creation from quark- antiquark annihila- 
tions, is given by cr ~ X^T^ /mf 2- The corresponding 
interaction rate exceeds the Hubble expansion rate pro- 
(4) vided that T ;> {X-'^mf^/Mpy/^. For A - 0(1) and 
7711,2 ~ O(TeV), this implies rapid equilibration of N at 
temperatures as low as T ^ itin ~ O(GeV). The comov- 
ing number density of TV decreases due to pair annihila- 
tion as the temperature drops below ttijv until thermal 
freeze-out of N annihilation. However, for 771 ^v ~ "Tip and 
70,1,2 ~ O(TeV), thermal freeze-out leads to overabun- 
dance of N according to the Lee- Weinberg bound ,1(^]. 

Therefore obtaining the correct relic abundance for N 
requires a non-thermal scenario of DM production. Late 
decay of a scalar S that reheats the universe to a low tem- 
peratures Tr <C O(GeV) can be the origin of non-thermal 
DM production. Interestingly, as we argued above, such 
a scenario is also required for successful baryogenesis in 
this model. In fact, the same decay processes that gen- 
erate baryon asymmetry also produce the DM candidate 
N. The relic density of DM particles thus produced is 
given by 
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N quanta, although produced when the temperature 
of the plasma is T^ <C tojv, are ultra relativistic since 
their initial energy is set by the mass of decaying particle 
ms- To determine whether N is cold or warm DM, one 
needs to see whether N can reach kinetic equilibrium 
with the plasma quickly enough. The interactions of N 
with M^, via exchange of ^1,2, also couple N to tt's. At 
temperatures m^r/S ^ Tj. <C O(GeV), tt's are abundant 
and kinetic equilibrium of N with the plasma is achieved 
through its interactions with tt's. At lower temperatures, 
T]. < TO,r/3, the number density of tt's is exponentially 
suppressed. In this case N interaction with the plasma 
occurs mainly through n exchange with the photons. 

We have checked that kinetic decoupling temperature 
of iV is ~ O(MeV), which is similar to the neutralino 



DM [17|, if mi,2, ms ~ ©(TeV) and A - 0(1). Achiev- 



ing kinetic equilibrium with the plasma very rapidly 
makes N a cold DM candidate. For much larger val- 
ues of TOi_2, "ms and/or much smaller values of A, kinetic 
decoupling temperature will be much higher. If it ex- 
ceeds Tr, kinetic equilibrium of N will not be achieved, 
in which case N can be warm DM. 



Explanation of the baryon-DM coincidence 
problem 



The interesting point to note from Eqs. ()6|7p is that 
the DM and baryon number densities have very similar 
functional dependence on the model parameters. As a 
result, they can come within the same order of magni- 
tude without making very special assumptions. For the 
natural choice of parameters |A| ^ |A'| ^ C(l), and CP- 
violating phases of 0(1), one indeed finds that nj^/s is 
larger than ijb by a factor of a few. Since m^r ~ rup, 
this is directly translated into similar relation between 
the DM and baryon energy densities, thus providing a 
natural explanation for the DM-baryon coincidence puz- 
zle. 

The model in Eq. (jlj therefore leads to a successful 
realization of the "Cladogenesis" scenario [lal- In this 
scenario the DM and baryon densities are mainly con- 
trolled by the dilution factor and branching ratios of the 
decaying field S. Assuming 0(1) values for A and A', one 
finds rjB and n^ / s in the ballpark of measured values 
due to the smallness of I5 and Bri^2- The exact numbers 
can then be obtained by minor adjustments of A and A' 
about their natural values. 

One can construct explicit models of a late-decaying 
scalar field S that has a small branching ratio to certain 
fields as discussed in [I3, |20| • In the case at hand, for 
example, one can assign lepton number —2 to S within 
a gauged U{1)l model. This allows S coupling to the 
RH neutrino Vfi via a lepton number conserving term 
hSvRVji. As usual, vr participates in the neutrino Dirac 
Yukawa coupling terms h^Hi/jiL. The U{1)l symmetry 
allows renormalizable couplings of the form 151^1X1^2^ 
between S and A'i_2, but not terms like SNN. Assuming 



that rn^^ < ms < Stti^^, the main decay mode of S is 
S —T' HvfiL. Here rn^^ is the Majorana mass of the RH 
neutrino, which is generated after spontaneous breaking 
of f/(l)L by some Higgs field. The S decay rate is sup- 
pressed ex {hhi^Y, in addition to the three-body phase 
space suppression. This can explain a small decay rate 
leading to a low reheat temperature Tj-, and hence a small 
value of Ys- Production of ^"1^2 occurs via 5-body final 
state decays S — > Xi^2X12HvrL, which result in very 
small values of Bri^2- 



IV. DIRECT DETECTION AND COLLIDER 
SIGNALS 

The DM candidate N interacts with nucleons through 
its coupling to the m^ mediated by X, see Eq. ([T]). How- 
ever, N only couples to a particular chirality of up- type 
quarks. As a result, there are no effective interactions 
of the form ('07v)'!/'Ar(V'g'0g)/"^x between the DM and 
quarks (V'tv and V'g are four-component fermions rep- 
resenting N and quark fields respectively). The spin- 
independent interactions with nucleons arise from twist- 
2 quark operators and one-loop diagrams that couple N 
to gluons [2l[ . The corresponding elastic scattering cross 
section will be suppressed ex m^ , which gives rise to 
CTsi < 10^1^ - 10^1^ pb for rnx -- O(TeV). This is con- 
siderably below the reach of upcoming experiments. 

The spin-dependent cross-section is only suppressed 
ex m^ as one has effective interactions of the form 
(V'Af757^V'Jv)(V'g757A«V'g)/™l:- This resuhs in ctsd <, 
10~^ — 10~^ pb, for mx ^ O(TeV), which is much below 
the bounds from current experiments [22|, as well as the 
upcoming detection experiments. It is also significantly 
below the current LHC bounds on ctsd [23], but may be 
within the LHC future reach |24| . 

Indirect detection signals from DM annihilation will 
also be negligible in this model. The reason being that 
because of the low DM mass, any signal will be over- 
whelmed by the astrophysical background. 

DM interactions with matter have a novel signature in 
this model that may be seen. The effective interaction 
Nu'^cTd'^ leads to baryon destroying inelastic scattering 
of N off nucleons, similar to the model in [25, which may 
have an appreciable rate for nucleon decay experiments. 

Moreover, this model can be observed at the LHC. 
The colored scalars Xi_2 can be pair produced and sub- 
sequently decay to a d'^cP pair or to Nu'^^ where N con- 
stitutes missing energy. The final state of X pair produc- 
tion can be 4 jets, 3 jets plus missing energy, 2 jets plus 
missing energy, etc. It is interesting to note that the final 
state with 4 jets does not involve any missing energy. 

Extending this model to a SUSY version is straightfor- 
ward and brings further interesting prospects. A minimal 
extension is done in the context of MSSM by introducing 
the SUSY partners of X scalars and N fermion, denoted 
by X and N respectively, plus two iso-singlet color-triplet 
superfields Ai^2 with hypercharge Y = —4/3 required 



by anomaly cancellation considerations. The Lagrangian 
£now in Eq. ([ij will be elevated to a superpotential Wncw 
to include all these fields and their interactions. 

If i?-parity is conserved, then the LSP, denoted by x, 
will be stable and can be a DM candidate as well. This 
leads to the possibility of multi-component DM consist- 
ing of TV and x- The LSP component of DM may be 
detected via direct searches even if it is subdominant. 
For example, x can be the Higgsino, or the superpart- 
ner of A'^ [2y] , both of which have elastic scattering cross 
sections within the reach of upcoming experiments. The 
multi-component DM scenario may also be detected via 
indirect searches. For example, a subdominant Higgsino 
component can lead to a detectable gamma-ray signal 
that is compatible with the Fermi-LAT bounds \0\. In 
this case, concordance between the collider, direct, and 
indirect signals would imply the necessity of a dominant 
non-LSP component of DM. However, complementary 
signals from direct and indirect searches are more con- 
strained due to the subdominance of the LSP compo- 
nent [23. 

In a SUSY set up one has to ensure that S decay does 
not lead to overproduction of the LSP component of DM. 
This can be achieved , for example, if S is an i?-parity 
even scalar that belongs to the visible sector [20]. We 
note that in i?-parity violating models N is the only DM 
candidate. In this case production of i?-parity odd parti- 
cles does not need to be suppressed since the LSP is not 
stable. 

At the LHC, the i?-parity odd colored fields X are pair 
produced and then decay to d'^d'^, if-N, u'^N final states. 
All of these final states involve jets plus missing energy. 
Also, if X is lighter than squarks, then the squark decay 
mode may contain the information of X. For example, 
u'^ — > AT + TV, X + N, with iV, X subsequently decaying 
into the LSP. The i?-parity even colored fields X decay 
to d'^d'^, Nu'^, d'^d'^, Nu'^ final states, with missing energy 
arising from A'' or from the LSP. An important difference 
is that the SUSY version of the model contains both 4 
jets and 4 jets plus missing energy final states. 

Finally, let us briefly mention some of the other phe- 
nomenological implications of this model. The Xd'^d'^ 
coupling can lead to K^-K^ or B^-B^^ mixing. However, 
color conservation does not allow any tree-level contribu- 
tion to the mixing term, and one-loop contributions that 
generate the relevant operators satisfy the experimen- 
tal constraints easily [Ij]. This model also gives rise to 
neutron- antineutron oscillations through the dimension-9 
operator G oc [u'^d'^s'^)'^ /{m\mN)- The oscillation time 
scale is given by t ~ 1/(2.510"^ G) s, where 2.5 x 10~^ is 
the value for the hadronic form factor [28|. The current 
bound t < 0.86 x 10^ sec ^ requires that G < 3 x lO'^s 
GeV~^. Since all possible flavor combinations appear in 
the expression for baryon asymmetry, see Eq. (jH]), it is 
easy to satisfy the oscillation bound while obtaining the 
correct baryon asymmetry 13 [. 



V. CONCLUSION 

In conclusion, we considered a simple extension of the 
SM that gives rise to baryogenesis and has a DM candi- 
date of ©(GeV) mass. Two colored scalars with O(TeV) 
mass and a singlet fermion are required in a minimal 
set up to generate the baryon asymmetry of the uni- 
verse via renormalizable baryon number violating inter- 
actions. The singlet fermion becomes stable and can play 
the role of a DM candidate, while avoiding rapid proton 
decay, when it is nearly degenerate in mass with the pro- 
ton. None of these explanations requires the existence of 
SUSY. 

DM and baryon asymmetry are produced non- 
thermally form out-of-equilibrium decay of new colored 
fields. Non-thermal baryogenesis is motivated when cou- 
plings associated with the new fields are not chosen to be 
artificially small. Obtaining the correct DM relic abun- 
dance for an ©(GeV) particle that interacts with matter 
via exchange of ©(TeV) scale fields also requires a non- 
thermal scenario. This mechanism correlates the DM 
abundance and baryon asymmetry for natural choice of 
parameters. Since DM mass is O(GeV), a correlation be- 
tween the number densities is directly translated into a 
similar relation between the energy densities. Therefore 
this model provides a natural explanation of the baryon- 
DM coincidence puzzle. 

The model predicts DM-nucleon scattering cross sec- 
tions crsi ^ 10-^*' - 10-1^ pb and ctsd <, 10"^ - 10"^ 
pb are far below the bounds from current and upcoming 
direct detection experiments. The predicted ctsd may 
however be within the LHC future reach. The indirect 
signals from DM annihilation are also negligible due to 
its low mass. However, the model may be probed at the 
LHC via the colored scalars. Pair production and subse- 
quent decay of these particles results in final states with 
N that constitutes missing energy. Final states with 4 
jets, 3 jet plus missing energy, and 2 jet plus missing en- 
ergy can arise but, interestingly, there will be no final 
states with 4 jets plus missing energy. 

Extending this model to a SUSY version is straight- 
forward. In i?-parity conserving models the LSP, which 
may be the scalar partner of iV, is also stable and can be 
a DM candidate too. This allows a scenario with multi- 
component DM, which has a better prospect for direct 
and indirect detection experiments. In the SUSY exten- 
sion, final states with 4 jets with missing energy can also 
arise. 
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